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The  construction  and  analysis  of  repair  models  is  an  important  area  in  reliability.  A 
commonly  used  model  is  the  minimal  repair  model.  Under  this  model,  repair  restores 
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models  that  could  be  classified  as  “better-than-minimal”.  Under  Kijima’s  models, 
the  system,  upon  repair,  is  functionally  the  same  as  a  working  system  of  lesser  age 
which  has  never  experienced  failure.  In  this  paper,  we  present  a  new  approach  to  the 
modelling  of  better-than-minimal  repair  models.  Using  this  approach,  we  construct  a 
general  repair  model  that  contains  Kijima’s  models  as  special  cases.  We  also  study  the 
problem  of  estimating  the  distribution  of  the  time  to  first  failure  of  a  system  maintained 
by  general  repair.  We  make  use  of  counting  processes  to  show  strong  consistency  of  the 
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1  Introduction 


Many  systems  are  maintained  and  kept  going  by  performing  repair  upon  each  failure.  The 
construction  and  analysis  of  plausible  repair  models  is  therefore  an  important  area  in  relia¬ 
bility.  Let  F  be  the  distribution  of  the  time  to  first  failure  of  the  system.  The  repair  models 
we  consider  in  this  paper  postulate  that  the  distribution  of  the  interfailure  times  depend 
in  some  way  on  F.  Early  works  on  repair  models  assume  that  repair  restores  the  state  of 
a  failed  system  to  a  level  equivalent  to  a  new  one  each  time.  This  is  the  so  called  perfect 
repair  model.  Clearly,  this  model  is  inadequate  to  model  most  repair  processes.  It  is  more 
reasonable  to  expect  the  distribution  of  the  remaining  life  to  vary  from  one  failure  time  to 
another.  One  such  model  often  used  in  the  literature  is  the  minimal  repair  model.  Under  this 
model,  repair  restores  the  system  to  its  state  just  before  failure.  Brown  and  Proschan  (1983) 
introduced  a  model  that  combines  both  perfect  and  minimal  repairs.  Under  their  model,  at 
the  time  of  each  repair  either  a  perfect  repair  occurs  with  probability  p  or  a  minimal  repair 
occurs  with  probability  1  —  p. 

There  has  been  a  search  for  repair  models  where  the  interfailure  times  are  stochaistically 
larger  than  in  the  case  of  the  minimal  repair  model.  Such  models  are  loosely  called  “better- 
than-minimal”  repair  models.  Kijima  (1989)  introduced  two  models  under  which  a  system, 
upon  repair,  is  functionally  the  same  as  an  identical  system  of  lesser  age.  Unlike  the  minimal 
repair  model,  Kijima’s  models  assume  that  the  distribution  of  interfailure  times  depend  on 
F  not  only  through  the  age  of  the  system  at  failure  but  on  the  degree  of  repair  as  well.  In 
this  paper  we  introduce  a  general  repair  model  in  which,  instead  of  reducing  the  effective 
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age  of  the  system  at  failure,  we  supplement  its  remaining  life  in  an  appropriate  way.  Such  a 
model  is  a  better-than-minimal  repair  model,  and  contains  Kijima’s  models  as  special  cases. 
We  will  also  study  the  problem  of  estimating  F  based  on  repair  data  on  n  identical  systems 
working  independently  and  maintained  using  this  general  repair  model. 

In  Section  2  we  describe  the  general  repair  model.  Section  3  discusses  stochastic  pro¬ 
cesses  needed  to  carry  out  the  estimation  process.  An  estimator  Fn  of  F  is  introduced  in 
Section  4  and  its  strong  consistency  is  proved.  The  estimator  of  F  is  the  solution  of  the 
integral  equation  Fn{')  —  io(l  —  Fn—')dAn  where  is  the  Aalen— Nelson  estimator  of  the 
hazard  function  A  of  F.  In  Section  5  the  weak  convergence  of  y/n{Fn  —  F)/F  to  a  Brownian 
Motion  is  established.  This  in  turn  is  used  to  construct  a  confidence  band  of  the  Hall-Wellner 
type  for  F„  in  Section  6. 

Estimation  of  the  distribution  of  the  time  to  first  failure  of  a  repairable  system  was  origi¬ 
nally  done  through  nonparametric  maximum  likelihood  techniques.  Whitaker  and  Samaniego 
(1989)  used  such  techniques  to  estimate  the  distribution  of  the  time  to  first  failure  of  a  system 
maintained  under  the  Brown-Proschan  model.  An  alternative  to  the  nonparametric  maxi¬ 
mum  likelihood  approach  is  the  use  of  the  theory  of  counting  processes.  Hollander,  Presnell, 
and  Sethuraman  (1992)  rederived  and  extended  the  large  sample  results  of  Whitaker  and 
Samaniego  using  such  an  approach,  making  use  of  the  methods  developed  for  the  analysis  of 
the  censored  life-data  model.  This  paper  extends  their  results  to  the  general  repair  model 
discussed  informally  in  this  introduction  and  formally  defined  in  Section  2. 
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2  Description  of  the  Model 


Prior  to  describing  the  model,  we  introduce  some  notation.  In  this  paper,  the  terms  increeis- 
ing  and  decreasing  will  be  used  loosely;  a  sequence  {aj}j>i  is  said  to  be  increasing  if  a,-  <  aj 
for  i  <  j  and  a  function  /(•)  is  said  to  be  increasing  if  f{x)  <  f{y)  for  x  <  y.  Similar 
definitions  apply  to  decreasing  sequences  and  functions.  When  we  require  the  condition  that 
a,-  <  Uj  for  i  <  j,  we  will  say  that  the  sequence  {aj}  is  strictly  increasing.  Similarly,  x  would 
be  called  positive  if  a:  >  0  and  negative  if  x  <  0.  The  notation  x  A  y  will  be  used  in  place 
of  min(x,i/).  Similarly,  x  V  y  will  be  used  in  place  of  max(x,y).  An  integral  of  the  form  /J 
would  denote  integration  on  (0,t].  We  use  the  Ito  integral  when  the  integrating  measure  is 
a  stochastic  process.  The  function  /(•)  will  denote  the  indicator  function,  so  that,  I{A)  =  1 
if  A  occurs  and  0  otherwise.  Finally,  we  will  adopt  the  convention  that  0/0  =  1. 

For  any  function  /,  we  will  use  f-{t)  for  lim^it  f{s)  and  f+{t)  for  limsjt  f{s).  A  function 
is  cadlag  if  it  is  right-continuous  and  has  left-hand  limits.  It  is  caglad  if  it  is  left-continuous 
with  right-hand  limits.  Let  T  <  oo.  We  will  denote  by  Z)[0,r]  the  cadlag  functions  on  [0,r] 
and  by  Z)_[0,r]  the  caglad  functions  on  [0,r].  Unless  otherwise  stated,  we  will  assume  the 
Skorohod  topology  for  Z)[0,  T].  The  supremum  norm  on  [0,  T]  will  be  denoted  by  ||  •  ||o  and 
“  =^”  will  denote  weak  convergence  of  probability  measures  on  D[0,  T]  under  the  Skorohod 
topology. 

/ 

Let  G  be  a  distribution  function.  We  will  denote  the  survival  function  1  —  G  by  G.  The 
hazard  function,  ff,  of  G  is  defined  by  H{-)  =  /q  dG/{l  —  G_).  The  distribution  G  is  said  to 


3 


be  IFR  if  G{t  +  x)IG{t)  is  decreasing  in  t  for  each  x.  Finally,  all  random  variables,  unless 
otherwise  stated,  will  be  assumed  to  be  defined  on  the  complete  probability  space  {Cl,^ 

We  will  now  formally  define  what  we  will  refer  to  as  the  general  repair  model.  For 
any  distribution  function  jF,  ^  6  (0,1]  and  a  €  [0,oo),  consider  the  family  of  distribution 
functions  F^{x)  =  F(6x  +  a)/F{a),  x  >  0  .  We  will  see  that  this  family  provides  us  with  a 
rich  class  of  distributions  for  the  remaining  life  of  a  given  system  subject  to  repair.  Consider  a 
system  put  into  operation  at  time  5*0  =  0  using  a  brand  new  unit  whose  life  distribution  is  F. 
Upon  each  failure,  the  system  is  repaired  in  negligible  time  and  put  back  into  operation.  Let 
denote  the  sequence  of  failure  times  of  the  system  and  let  Tj  =  Sj  —  Sj-i,  i  >  1,  be 
the  interfailure  times.  We  further  assume  that  the  Tj’s  are  strictly  positive.  A  repair  model 
describes  the  joint  distribution  of  the  random  variables  {Tj}.  In  this  paper,  we  describe  a 
general  repair  model,  based  on  two  sequences  and  called  the  effective  ages 

and  life  supplements,  respectively,  satisfying 

=  0, 01  =  1,  Aj  >  0,  6j  €  (0, 1]  and  Aj  <  Aj.i  +  ej-iTj.i  for  i  >  1  .  (2.1) 

The  model  is  obtained  by  specifying  the  joint  distributions  of  the  interfailure  times  {Tj}  as 
follows: 

P{Ti<t\Ai,ei,T, . =  hit>0,j>l  .  (2.2) 

The  family  of  distributions  }  are  stochastically  ordered  in  6,  i.e.,  0  <  6  implies  that 
F^  >  F^  ,  for  each  a  (i.e.,  F^(t)  <  F^  (t)  for  every  t).  Thus  for  the  general  repair  model 
described  by  (2.1)  and  (2.2),  the  distribution  of  Tj  which  is  Ff^  is  stochastically  larger 
than  F\.,  i.e.,  is  better  than  a  working  system  of  age  Aj.  Furthermore,  from  (2.1)  we  can 
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see  that  for  each  i  >  1  the  effective  age,  Aj+i ,  of  the  system  after  the  jth.  repair,  is  less 
than  its  effective  age  ,  Aj  +  QjTj,  just  before  the  jth  failure  which  in  turn  is  less  than  the 
actual  age  Sj.  Thus  the  general  repair  model  defined  by  (2.1)  and  (2.2)  can  be  considered 
as  a  better-than-minimal  repair  model  and  as  we  shall  see,  contains  the  perfect  repair  and 
minimal  repair  models.  We  will  illustrate  this  general  repair  model  and  the  terms  effective 
ages  and  life  supplements  through  examples  discussed  in  the  succeeding  paragraphs. 

Consider  the  case  when  0^  =  1  and  Aj  =  0  for  j  >  1.  Then  (2.1)  is  automatically 
satisfied  and  (2.2)  reduces  to 

P{Tj  <  ilTi, . . . ,  r,_i)  =  Fj(t)  =  F{t).  (2.3) 

From  this  we  see  that  the  Tj’s  are  independent  with  common  distribution  F.  This  corre¬ 
sponds  to  the  perfect  repair  model. 

Next,  consider  the  case  when  0j  =  1  and  Aj  =  Sj-i  for  each  j.  Clearly,  Aj+i  =  Sj  — 
Aj  -|-  QjTj  and,  hence,  (2.1)  is  satisfied.  Moreover,  under  these  conditions,  (2.2)  reduces  to 

=  (2.4) 

Hence,  we  see  that  this  case  corresponds  to  the  minimal  repair  model. 

We  will  now  show  that  Kijima’s  models  can  be  derived  through  suitable  choices  of 
{Aj}  and  {0j}  satisfying  (2.1)  and  (2.2).  Leb{Dj}j>i  be  a  sequence  of  random  variables 
independently  distributed  on  [0, 1]  and  independent  of  other  processes.  Consider  the  case 
when  Qj  =  1  for  each  j  and  Aj  =  DiTi  for  j  >  1.  Since  Ajj^i  —  Aj  +  DjTj  and 
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Dj  <l  =  0j  for  each  j,  then  (2.1)  is  satisfied  and  (2.2)  reduces  to 

p{Ti  <  i|r„  A,  I  <  i  <  i  - 1)  =  ■ 


(2.5) 


This  is  Kijima’s  Model  I.  In  this  model,  upon  the  {j  -  l)th  repair,  the  time  to  next  failure, 
Tj,  of  the  system  has  the  same  distribution  as  the  time  to  first  failure  of  a  system  whose  life 
distribution  is  F\..  Hence,  upon  repair,  the  system  whose  actual  age  is  Sj  is  functionally  the 
same  as  an  identical  system  of  age  Aj+i  which  has  never  experienced  failure.  This  explains 
the  use  of  the  term  effective  age. 

Consider  the  case  when  0j  =  1  for  each  j  and  Aj  —  A)2fc  for  ;  >  1.  Since 

Aj+i  =  Dj{Aj  +  Tj)  and  Dj  <  1  =  Qj  for  each  j,  then  (2.1)  is  satisfied  and  (2.2)  reduces 
to  (2.5).  This  is  Kijima’s  Model  11.  Moreover,  when  Dj  is  1  with  probability  p  and  0  with 
probability  1  —  p,  we  obtain  the  Brown- Proschan  model. 

Up  to  this  point  we  have  restricted  the  0j’s  to  be  identically  equal  to  1.  We  will 
now  describe  repair  models  obtained  through  other  choices  of  the  life  supplement  sequence 
{0j}.  Recall  that  Ff  >  if  6  <  6' .  This  implies  that  the  smaller  6  is  the  larger  is  the 
expected  remaining  life  of  the  system.  Hence,  we  can  use  ^  as  a  measure  of  how  repair 
supplements  the  expected  remaining  life  of  the  system.  This  explains  the  use  of  the  term 
life  supplement.  If  a  minimal  repair  was  performed  at  the  time  of  the  first  failure,  then  T2 
would  have  the  distribution  .  If  we  want  a  longer  expected  life  for  then  we  can  use 
the  distribution  F^^  for  some  02  satisfying  0  <  02  <  1-  Starting  with  the  distribution  F^^ 
for  T2  and  using  minimal  repair  upon  the  second  failure,  the  random  variable  T3  would  have 
the  distribution  F\^  where  A3  =  Ti  +  Q2T2.  Again,  if  we  want  a  longer  expected  life  for 
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Ts  we  can  use  the  distribution  F®/  for  some  63  satisfying  0  <  ©3  <  1.  With  this  in  mind, 
we  define  ©i  =  1,  Aj  =  ZCi  ^iTi  and  0  <  ©j  <  1  for  j  >  1.  Under  these  conditions, 
Aj+i  —  Aj  +  QjTj  for  j  >  1,  satisfying  (2.1),  and  (2.2)  reduces  to 

P(Ti  <  (ITi,  ©i,  1  <  i  <  J  -  1)  =  F^‘(t)  .  (2.6) 

We  will  refer  to  this  as  the  “supplemented-life”  repair  model.  By  definition,  the  system 
enjoys  a  larger  expected  remaining  life  under  a  supplemented-life  repair  than  it  would  under 
a  minimal  repair. 


3  Some  Fundamental  Processes 

There  is  a  close  connection  between  repair  models  and  censored  life-data  models  that  enables 
one  to  use  techniques  developed  for  censored  life-data  models  in  the  analysis  of  repair  models. 
Fix  a  T  >  0  and  define  the  processes  N  and  Y  by 
N{t)  =  sup{j  >  1  :  Yj-  <  i,  Sj  <  r} 

and 

Y{t)  =  sup{i  >l:Aj<t<  {Xj  A  [Aj  +  Qj{T  -  S'j-i)])} 
where  Xj  =  Aj  +  OjTj.  Let  6j  =  I{Sj  <  T)  and  Xj  =  Xj  A  [Aj  +  ©j(r  -  5,_i)].  Then 
the  random  variables  {(Yi,^i),  (X2,  (^2), . . .}  can  be  viewed  as  observations  coming  from  a 
censored  life-data  model.  A  general  repair  model  observed  during  a  period  of  length  T  is 
akin  to  a  survival  study  where  a  subject  j  enters  the  study  at  age  Aj  and  either  dies  during 
the  study  at  age  Xj  or  leaves  the  study  by  age  Aj  -|-  Qj{T  —  Sj-i).  So  that  at  the  completion 
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of  the  study  only  the  variables  {(Xj,  ^i),  (A'2,  ^2),  •  •  •}  are  actually  observed.  The  random 
variable  N(t)  represents  the  number  of  observed  (uncensored)  deaths  by  time  t  and  Y{t) 
the  size  of  the  risk  set  at  time  t.  Since  the  processes  N  and  Y  are  fundamental  in  the 
estimation  of  the  survival  function  in  a  censored  life-data  model,  it  is  reasonable  to  expect 
these  processes  to  play  a  similar  role  in  estimating  F  under  general  repair. 

Let  A  be  the  hazard  function  of  F  and  define  the  process  M  =  N  —  f  Y dA.  The  process 
M  plays  an  important  role  in  establishing  the  large  sample  properties  of  the  Whitaker— 
Samaniego  estimator.  Based  on  this  observation,  it  is  reasonable  to  expect  M  to  play  a 
similar  role  in  establishing  large  sample  properties  of  Fn-  We  now  find  expressions  for  the 
mean  and  covariance  functions  of  M.  To  simplify  expressions  we  shall  use  the  following 
notation.  For  j  >  1,  let 


Nj{t)  =  /(5,<t), 

ffj(s)  =  I{Aj  +  Qj{s-Sj.r)<t)  , 

Gj{t)  =  I{Sj-^<t<Sj)  , 

=  A{Aj  +  Qj{t  —  Sj-i))  , 

Klis)  =  /(4  <s<A,  +  0,(f-5,_i))  , 


and 

Mj{t)  =  Njit)-f^G^dAj  . 

We  can  now  rewrite  M  as 


=  Y.r  H‘{3)dNi(s) 


5,-; 


H]{s)Gi{s)dA,{s) 
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=  E‘iXi<t,Si<T]- 

j 

HA,  <  3  <  Xi  A  [Ai  +  0,(r  -  S,-i)J)rfA(i) 

=  x:  Njimasi)  - 1;  [  Kr^‘(3)dA(s) 

=  N{t)  -  rY{s)dA{s)  =  M(t).  (3.1) 

Jo 

Hence  to  evaluate  the  mean  and  covariance  functions  of  M  it  is  enough  to  evaluate  those  of 
Sl_,  Hi{3)dMi(3). 


Lemma  3.1  For  fixed  t  and  t' , 


0 

for  all  j 

(3.2) 

0 

for  i  ^  j 

(3.3) 

E[ 

Jo 

for  all  j. 

(3.4) 

Proof:  For  each  j,  let  Fj  =  (7{Aj,Qj-,Ti, . . .  The  idea  behind  the  proof  is  quite 

simple.  Upon  conditioning  on  F’j,  we  are  left  with  expressions  of  the  form 

J g{s-,Aj,Qj,Sj-i)P{Xj  >  h{s-,  Aj,Qj,  Sj-i)\Xj)dA{s) 

for  some  functions  g,h  .  These  expressions  are  then  easily  simplified  using  (2.2).  We  now 
carry  out  the  calculations. 

To  prove  (3.2),  note  that 

Elfi  H'GjdAilX,]  =  £{[' Kj''^’dA\r,] 
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=  J‘Kf{s)P{Xi>s\ri)dA{s) 

I  ’  ^  h(Ai) 

F(iA|/i,  +  6j(r-s,-i)] 

f{Ai) 


El  HidNilEj] 

J  5j_i 


This  proves  (3.2).  From  this  it  follows  that  for  i  <  j 


E[r  HldMiT  H}dMj\Ej]=  r  HldMi-E[r  H}dMj\Ej]^0. 


This  proves  (3.3).  To  prove  (3.4),  note  that 


Now, 


;|f  H'diii  r  H'dMiiTj] 

=  E[Ni{T)H‘{Sj)Hf  (Si)\Ei] 

T 

-  ElNi(T)Hl(Si) I  HiG,dAi\ri] 


ElNi{T)H‘  (Si 


H)G,dAj\Ei] 


-  ElJ^  B‘jGjdAj  H' GjdAilEj] 


ElNj(T}H‘{Sj)H'  (5,)|r,l 
=  P{X,  <  (i  A  (')  A  lA,  +  ej{T  -  s,-,)i|?-,) 

,  ^[(iA(')A(4,  +  e,(r- v.))l 


F{Ai) 


=  r 
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and 


E[Ni(T)H'j{Si)  £  HiGidAi\ri] 

t 

=  Kj{s)I{s  <Xj<tA  [Aj  +  0j(T  -  Sj.i)])dA{s)\Xj] 

rtAt 

=  jf  Kj{s)P(s<Xi<tAlAi  +  ej{T-Si-iW,)dA(s) 

=  L  T(Ai) 


By  symmetry, 


Finally, 


E[Nj(T)Hi  (5,)  £  H'GidAi \T,] 

=  r'  ^  -  ^‘-'"^dAis). 

Jo  ^ 


Ei£  H'GidAj  £  H'-  GjdAj  I 
=  E[j'  [  Kf  {s}Kf  (u)I(Xj  >sV  u)dA(s)dA(u)\rj] 
—  f  f  Kf  {s)Kj{u)P[Xj  >  s  \/ u\Xj]dA{s)dA{u) 


+ 


rtAt 

L 


F.(s) 


AA(s)<JA(s) 


■F(Ai) 

„T,  ,nt'  A  lAi  +  e,(r-5,-,)])-.F(u) 

F{Ai) 


rtAt 
=  /„ 


dA(u) 


+ 


F((  A  lAi  +  0,(r  -  5,-i)])  -  F(s) 


F(Aj) 


rtAt 

/„ 


dA(5) 


(3.7) 


(3.8) 


(3.9) 
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Substituting  (3.6),  (3.7),  (3.8),  and  (3.9)  into  (3.5)  we  get 

EiT  HjdMi  r  H‘dMi\Ei] 

rt/\t  _ 

=  /  /fJ(s)P[Jif,  >s|:7)|[l-AA(s)]<iA(s) 

*/ 0 

/ 

=  BIjf'"  /<r*'(s)(l  -  AA(s))<iA(s)|^)) 

This  proves  (3.4)  and,  thus,  completing  the  proof.  <C> 

From  this  result  we  obtain  the  following  expressions  for  the  mean  and  covariance  func¬ 
tions  of  M. 


Theorem  3.1 


EM 

o 

II 

< 

1 

II 

(3.10) 

Cov{M{t)M{t)) 

rtAt 

=  /  EY{1  -  AA)dA 

Jo 

(3.11) 

Proof:  From  (3.1)  and  (3.2)  it  follows  that 

EM  =  Ej2r  H^j{s)dMj{s)  =  0. 
j 

This  proves  (3.10).  To  prove  (3.11)  we  use  equations  (3.1),  (3.3),  and  (3.4)  to  obtain  the 
following 


T  jP  / 

E{M{t)M(t'))  =  EJ2  I  HldMi  [  nidMj 
=  Kj''"^il-AA)dA 

i-tAt' 

=  /  EY{1  -  AA)dA. 

Jo 
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This  completes  the  proof.  <C> 


4  The  Estimator  and  Its  Uniform  Consistency 

Suppose  that  we  observe  n  independent  copies  of  the  processes  N  and  F  on  a  finite  interval 
[0,  T].  We  will  look  at  the  problem  of  estimating  F  based  on  these  observations.  Throughout 
the  rest  of  the  paper,  we  will  assume  that  (i)F(T)  <  1  and  (ii)  F  is  IFR  .  Assumption  (i) 
allows  us  to  bound  EY(t)  away  from  zero  uniformly  on  [0,3^].  This  is  pivotal  in  proving 
uniform  consistency  of  the  estimator.  Define  N{t)  =  sup{j  >  1  Sj  <  t}.  The  random 
variable  N{t)  is  the  number  of  system  failures  by  time  t.  Let  N*{t)  —  sup{7  >  1  :  S’-  <  t} 
where  {5^  }  is  the  minimal  repair  process  based  on  F.  Assumption  (ii)  gives  us 

p{si  -  st_,  >  xisr_,  =  ()  > 

for  all  jfc  >  1  and  t  €  [0,  T].  It  follows  that  the  evaluation  at  time  t  of  a  renewal  process  having 
recurrence  time  distribution  G{x)  =  1  —  F{T  +  x)fF{T)  stochastically  dominates  N*{t)  on 
[0,  T]  which  in  turn  stochastically  dominates  N{t).  (The  counting  process  Ni  stochastically 
dominates  N2  on  [0,r]  if  P(Ni{t)  >  n)  >  P{N2{t)  >  n)  for  all  n  and  t  €  [0, T].)  This 
guarantees  that  N{T)  has  finite  moments  of  all  orders.  We  will  need  this  in  proving  weak 
convergence. 

Let  Nn  and  Yn  be  the  sum  of  the  first  n  independent  copies  of  N  and  Y,  respectively. 
We  note  that  from  the  viewpoint  of  a  life-data  model,  dF{t)f(l  —  F^{t))  represents  the 
instantaneous  failure  rate  at  time  t.  A  straightforward  way  of  estimating  this  failure  rate  is 
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by  using  the  ratio  of  the  number  of  observed  deaths  at  time  t  to  the  number  at  risk  at  time 
t.  Hence,  a  natural  estimator  of  the  hazard  function  A  is  given  by 


Lit)  =  / 

Jo 


*  Jn  •  dNn 
Yn 


where  Jn(t)  =  /(K(t)  >  0)  for  t  €  (0,r].  The  estimator  A„  is  referred  to  as  the  Aalen- 
Nelson  estimator  of  A.  Since  F{t)  =  /o(l  -  F-)dA  ,  it  is  reasonable  to  require  an  estimator 
F„  of  F  to  satisfy 

Fnit)=  r(l-A-)dA„.  (4.1) 

Equation  (4.1)  is  sometimes  referred  to  as  a  Volterra  integral  equation.  Its  solution  is  given 

by 

A(t)  =  11(1  -  dLis))  =  11(1  -  AAn(5)) 

$<t  s<t 

where  n5<t(l  dAn{s))  denotes  the  product  integral  (see  Gill  and  Johansen  (1990)). 


The  rest  of  the  paper  would  now  be  devoted  to  a  study  of  the  large  sample  properties 
of  jPn*  Under  a  minimal  repair  assumption,  the  process  M  of  Section  3  turns  out  to  be 
a  martingale.  This  enabled  Hollander,  Presnell,  and  Sethuraman  (1992)  to  use  martingale 
techniques  in  proving  large  sample  results  for  the  Whitaker—Samaniego  estimator.  For  the 
general  repair  model,  finding  a  suitable  filtration  with  respect  to  which  M  is  a  martingale 
proves  to  be  a  formidable  task.  A  lot  of  the  difficulty  is  due  to  the  fact  that,  in  general, 
the  Xj'^s  are  not  monotonic.  This  makes  the  existing  methods  for  computing  a  compensator 

if. 

inapplicable  in  the  general  case.  One  could  try  working  with  the  failure  process,  say  N, 
associated  with  the  Sj^s  instead  of  the  failure  process  N  associated  with  the  Xj^s.  Since 
the  Sj's  are  strictly  increasing  it  would  not  be  too  hard  to  find  the  compensator  of  N.  The 
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problem  would  then  turn  to  transforming  results  on  N  into  results  on  N.  Unfortunately,  we 
found  this  to  be  an  equally  formidable  task.  However,  a  closer  look  at  Theorem  3.1  reveals 
that  although  M  may  not  be  a  martingale  with  respect  to  the  history  of  N,  it,  nevertheless, 
exhibits  the  same  mean  and  covariance  structure  it  would  have  if  it  were  a  martingale. 
Fortunately  these  features,  in  conjunction  with  techniques  used  by  Gill  (1980)  for  Markov 
renewal  processes,  are  sufficient  to  prove  large  sample  results. 

Since  most  of  the  distributions  considered  in  reliability  are  continuous,  we  will  assume 
throughout  the  rest  of  the  paper  that  F  is  continuous.  This  will  help  us  avoid  unnecessary 
complications.  This  by  no  means  limit  our  results  to  the  continuous  case.  Most  of  our 
arguments  carry  over  to  the  discontinuous  case  with  very  little,  if  any,  modifications. 

The  process  Y  can  be  written  as  the  difference  of  two  left-continuous  increasing  pro¬ 
cesses  Yi ,  Y2  defined  by 

ri(i)  =  sup{j  >l:Aj<  <  T} 

and 

Y^it)  =  sup{i  >  1  :  X,-  <  f,  <  T} 

where  Xj  -  Xj  A  [Aj  -f  Qj{T  -  Let  Fin  and  Y2n  be  the  sum  of  the  first  n  independent 

realizations  of  Yi  and  Y2,  respectively.  Then  the  process  Fn  can  be  considered  as  the  result 
of  the  three  mappings 

^  Flu  ^ 
n  n  '  n 

going  through  the  spaces 

£»[0,r]  X  D.[0,T]  X  D_[0,r]  — >  D[0,r]  X  D_[0,r]  — ^  D[0,T]  D[0,T]. 


Yn 


An  - >  Fn 
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Under  the  supremum  norm  H  •  ||,  the  first  mapping  is  clearly  continuous.  To  show  that  the 
second  mapping  is  continuous  we  use  the  following  result. 


Lemma  4.1  Let  =  Nnln.Y^  =  Yn/n  and  a  =  inf(o,T]  EY{t)  .  Then  with  probability  1 


An  j 


dEN 


EY 


< 


iK-EYW^-jENm  +  WK-ENg] 
a{a  -  UK  -  EY\\^) 

+hK-ENg 

a 


(4.3) 


for  sufficiently  large  n. 

Proof:  Since  F{T)  <  1  then  P{S\  >  T)  >  0.  It  follows  that  for  t  €  (0,  T] 
EYit)  >  P{Y{t)  >  0)  >  P{Y{T)  >  0)  >  P(5i  >  T)  >  0. 


Hence,  a  >  0  and 


This  implies  that 


/ 


dEN 


EY 


< 


< 


i|lK  -  EYg  ■ 

a 


f  rf[(jv;  -  EN)  +  £iV] 
/(0,.l  ey-(ey-y:) 


+  |iiiv:  - 

IlK  -  £y|l?  ■  1  +  IK  -  ■EJvig  I 

cc(c  -  UK  -  EY\\l) 


+  IwK 


EN\ 


0 


T 

0 
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Now,  by  the  Glivenko-Cantelli  theorem 


max  • 


N 

^-EN 

n 


^-EYi 

n 


T  ny 


2n 

n 


-EY2 


0. 


From  (4.3)  and  the  continuity  of  the  first  mapping  in  (4.2),  we  get 


A 


dEN 


EY 


0  w.p.  1  . 


But  from  (3.10)  with  t  €  (0,r] 


I 


t  dEN 
EY 


A{t). 


Hence  , 


||An  —  A||o  — >  0  w.p.  1  as  n  — 00 


(4.4) 


Uniform  consistency  now  follows  from  (4.4)  and  the  continuity  of  the  product  integral  as 
shown  in  the  next  result. 


Theorem  4.1  \\Fn{t)  -  F(i)||^  — >  0  w.p.  1  as  n  — >  00. 

Proof:  Since  A„(t)  is  increasing  in  t  for  each  n  and  A„(T)  — >  it  follows  that 

lim  sup  |A„|((0,  rj)  <  00  w.p.  1.  This  guarantees  the  continuity,  under  the  supremum  norm, 
of  the  last  mapping  on  (4.2).  (See  Theorem  7  of  Gill  and  Johansen  (1990).)  <) 


5  Results  on  Weak  Convergence 


Define  =  \/n[-j^n  —  E]/F  and  the  process  M„  by 

Mn  =  Nn-J  F„dA. 
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Now,  by  Lemma  7.2.1  of  Shorack  and  Wellner 


hit)  -  F{t) 

m 


-  A- 

l-F 


F 


•d[A„-A] 


for  all  t  6  (Ojrj.  Hence 


where  Wn  =  Note  that  from  the  CLT,  the  finite-dimensional  distributions  of  Wn 

converge  to  that  of  a  Gaussian  process.  This  suggests  that  it  might  be  possible  to  obtain  a 
weak  convergence  result  for  Wn  and,  consequently,  for  Zn  in  view  of  (5.1). 


Theorem  5,1  Let  W  be  a  zero  mean  Gaussian  process  with  independent  increments  and 
variance  function  Var{W{t))  =  /q  EYdA,  Then  Wn  =>  W  on  D[0^T]  . 

Proof:  This  theorem  is  analogous  to  that  of  Lemma  3  in  Gill  (1980)  and  is  proved  sim¬ 
ilarly.  Note  that  Mn  is  the  sum  of  n  independent  copies  of  the  process  M.  That  the 
finite-dimensional  distributions  of  Wn  converge  to  that  of  W  follows  from  the  CLT.  It  re¬ 
mains  to  prove  tightness.  To  prove  this  it  suffices  to  show  that  there  exist  a  nondecreasing 
right-continuous  function  G  and  rj  >  1/2  such  that 

n-^E[iMM  -  MnimMr^it)  -  M„(.))^]  <  (G(«)  -  GimG{t)  -  Gis))^ 

foT0<s<t<u<T  (see  Billingsley  (1968)  page  133).  For  notational  convenience  denote 
A,tH  =  H{t)  -  His).  Then 

n-^EiAt^Mr,-A,tMnf  =  n-'E(Ai„M  • 
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n 

^2-'^[E{AtuM-AstM)]'^. 

n 


(5.2) 


Recall  that  N{t)  =  sup{j  >  1  :  5j  <  t}.  Let  V  =  iV  +  1  and  J5  >  (1  -  F{T))  L  Then 

EiAtM-AstM)  =  0  (5.3) 

E{AtM?  = 

<  B  ■  EV{T)  •  {At^Ff  (5.4) 

for  any  t  <u  and  0  <  a  <  1. 

\AtuM-A,tM\  <  {AtuN  •A,tN)  +  {l"YdA){jydA) 

-f-  •  J  Y dA  AgtN  •  J  Y dA 

<  1(1  +  B)V{T)Y  (5.5) 

From  (5.3)  and  Holder’s  inequality  it  follows  that  for  Or  =  1  —  €  (0, 1) 

E{AtuN  •  AstN)  +  E{£  YdA  •  J'  YdA) 

=  E{At^N- jydA)i-E{AstN- J\dA) 

<  B-EiAt^N-V{T)){AstF) 

+  B  •  E{AstN  ■  V(T))(AtuF) 

<  B  ■  [EV{Ty+^/^^Y  [E{AtuN)y^  (AstFy^ 

+  B  ■  [EV{TY+'^^^^yyE{A,tN)]^yAtuFyK 
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Now  E{AtuN)  <  B  •  {EV{T)){AtuF)  for  t  <  u.  Hence, 

E{\AtuM  ■  A,,M|)  <  Cx(At„F  •  A.^F^  (5-6) 

for  some  constant  Ci  provided  EV{Tf'^^l^i^  <  oo.  From  (5.5),  (5.6),  and  Holder’s  inequality 
it  follows  that  for  some  0:2  =  1  —  /?2  €  (0, 1)  and  constant  C2, 

E{AtuM-A,tMf  =  E{\AtuM  ■  A,tM)\^^\At^M  ■  AstM)\^'^^) 

<  (1  +  •  (F|At„M  •  AstM\Y^{EV{Tf^'^<^^')<^ 

<  C2(At„F  •  (5.7) 

provided  EViTY'^'^^^^  <  00.  Since  N(T)  has  finite  moments  we  can  choose  /?i  <  1/6  and 
/?2  <  2/5  so  that  aiQ:2  >  1/2.  Substituting  (5.3),  (5.4),  and  (5.7)  into  (5.2),  we  get 

n-^EiAtuMn  ■  AstMnf  <  C  ■  (A<„F  •  A.^F)’' 

for  some  constant  C  and  77  >  1/2.  This  completes  the  proof.  <) 

From  the  Skorohod-Dudley-Wichura  theorem  (see,  e.g.,  Shorack  and  Wellner  (1986) 
Theorem  2.3.4)  there  exist  processes  W'^  and  W  having  the  same  distribution  as  Wn  and 
respectively,  such  that 

psiK^'iv')  0 

where  ps  is  the  Skorohod  metric  on  F[0,  T].  In  light  of  this,  we  will  assume  that  we  actually 
have  ps{Wn,  W)  — >■  0.  To  prove  convergence  results  for  J\r„,  however,  we  would  be  needing 
uniform  convergence.  If  we  can  show  that  W  has  almost  surely  continuous  paths  then  we 
would  have  uniform  convergence  since  under  such  a  condition  the  supremum  norm  distance 
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is  equivalent  to  Skorohod  distance.  That  W  has  almost  surely  continuous  paths  is  shown  in 
the  next  theorem. 

Theorem  5.2  ||W„(t)  -  W{t)\\‘^  — ^  0  w.p.  1  as  n  — ^  oo. 

Proof:  From  Theorem  3.1  and  Chebyshev’s  inequality,  we  get 

P{\M{t)  -  M{to)\  >e)<  I  0  as  |t  -  to|  ^  0 

for  each  e  >  0  and  to  G  [0,  T],  i.e.,  M  is  stochastically  continuous  on  [0,  T].  This  implies  that 
W  has  almost  surely  continuous  sample  paths  (see  Theorem  2  of  Hahn  (1978)).  Hence  the 
result.  <C> 

Remark:  If  F  is  not  continuous  then  M  may  not  be  stochastically  continuous.  Hence, 
Theorem  5.2  may  fail  to  hold.  Fortunately,  we  can  construct  versions  of  Wn  and  W  so  that 
Theorem  5.2  holds  for  these  versions.  (See  Lemma  4  in  Gill  (1980).) 

Theorem  5.3  Define  Z{t)  =  jQ{EY)~^dW .  Then  ||Z„  —  Z\\o  — >  0  w.p.  1. 

Proof:  Let  =  {Fn-/F)  ■  [nJ„/y„]“L  It  follows  from  (5.1)  that  for  t  €  (0,r],  Z„(f)  = 
flHndWn.  The  result  now  follows  since  ||ifn  —  (jBF)“^  j|o  — >  0  w.p.  1  and  Hl'Fn  — IF||o  — ^ 
0  w.p.  1  (see  Lemma  5  of  Gill  (1980))  .  ^ 

Now  since  the  supremum  norm  distance  is  larger  than  the  Skorohod  distance,  ps,  it 
follows  that  ps{Zn,  Z)  — >•  0  w.p.  1  .  This  gives  us  the  following  result. 

Theorem  5.4  =  ^/n{Fn  —  F)fF  Z  on  D[0,T]  . 


21 


Proof:  See  Corollary  2.3.1  of  Shorack  and  Wellner  (1986)  .  0 


Corollary  5.1  Let  B  denote  Brownian  motion  on  [0,  oo)  and 


^W  =  /o' 


1 


EY  1 


dF 

-F  ■ 


Then  y/n{Fn  —  F)fF  B{C)  on  D[0,T]. 


Proof:  Note  that  {W^(<)}t€[o,r]  is  a  square  integrable  martingale  with  respect  to  Ft 
:  s  <t).  Furthermore,  it  is  easily  checked  that  {W)t  =  Jq  EYdA.  Hence 

Cov(Z(s),Z(t))  =  E[Z(sAt)f 

=  E  r\EY)-H{W) 

Jo 

fsAt 

=  /  (EY)-'^dA  =  C{s  At)  . 

Jo 

This  proves  the  result.  0 

Corollary  5.2  Let  B°  denote  a  Brownian  bridge  on  [0, 1]  and  K  =  C/{1  +  C).  Then 

V^j{Fn-F)=:^B%K)onD[0,T]  . 

Proof:  Let 

B°{t)  =  (1  -  i)  •  b{y^)  for  t  e  [o,/<:(r)]. 

The  result  follows  immediately. 
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6  A  Hall-Wellner  Type  Confidence  Band 


In  this  section  we  construct  a  confidence  band  for  F  similar  to  the  bands  of  Hall  and 
Wellner  (1980)  for  the  censored  life-data  model.  We  then  present  results  of  some  simulation 
studies  of  the  coverage  probabilities  of  the  band  using  repair  models  discussed  in  Section  2. 


For  t  in  [0,  T],  define  Ln{t)  =  I{Fn{t)  <  1)  and  let 

n  (f  \  —  ^nLn  dFn 
^  Jo  Ynin  '  l-Fn 


Define 


Kn{t)  = 


Cn{t) 

1  +  Cn{t) 


For  t  such  that  F„(<)  =  1  set  Kn{t)  =  1-  Then  the  result  of  Corollary  5.2  suggests  a 
confidence  band  for  F  of  the  form 


(6.1) 


where  is  such  that  P{suptg[o,i]  ^  =  1  -  a.  This  gives  our  band  an  asymptotic 

confidence  level  of  at  least  100(1  —  a)%.  To  justify  this  band  we  need  the  following  result. 


Lemma  6.1  ||/^„  —  0  w.p.  1  . 


Proof:  For  t  €  [OjT’j, 


\Cn{t)-C{t)\  < 


Jo  [YJn  .  (1  ■ 


(1_A)  EY-{1-F)\ 


dFn 
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Note  that  a  =  inftg(o,r)  >  0  (Lemma  5.1)  and  that  suptg(o,r)  ■^(^)  ^  These 

together  with  the  uniform  consistency  of  Yn/n  and  Fn  implies  that  JnL„  =  1  on  (0,r]  with 
probability  1  for  sufficiently  large  n.  Hence,  with  probability  1 

\C4t)-C{t)\  <  •  rn/"  •  (1  - -ft.)  - -ey  •  (1  - -F’)ll? 

1 

EY-{1-  F) 

for  sufficiently  large  n.  Since  Yn/n  and  Fn  are  uniformly  consistent  then  so  is  Cn-  The  result 
follows.  <C> 

This  leads  to  the  next  theorem  which  justifies  our  confidence  band. 

Theorem  6.1  A/n^(F„  —  F)  B°{K)  on  jD[0,r]  . 

Fn 

Proof:  The  result  follows  from  the  identity 

=  •  (v5(A  -  F)] .  <> 

F  ^  F.^  ^ 

Remarks: 

(i)  In  practice,  it  may  be  the  case  that  the  data  obtained  leads  to  Fn{to)  =  1  for  some 
0  <  to  <  T.  When  this  happens,  the  data  obtained  gives  us  a  confidence  band  only  on  the 
interval  [0, (t)  where  a  =  inf{t  G  [0, T]  :  Fn{t)  =  1}. 

(ii)  Let  X(i),  ^(2), . . . , -^(r)  be  the  distinct  ordered  values  of  the  X’s  whose  corresponding 
failure  times  are  within  [Oj^r].  Also,  let  6j  be  the  number  of  observations  with  value 
Then  for  computational  purposes  we  note  that 
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h{t)=  n  1 


YA^U)). 


and 


f)  (i\  —  „  ^ni^jj-l)) 

Dn{t)  —  n  2^  p.  . 


We  illustrate  the  band  using  examples  from  Section  2  and  the  gamma  and  Weibull 
distributions.  For  our  first  repair  model  we  will  use  Kijima’s  Model  II  where  the  D/s  are 
taken  to  be  uniformly  distributed  on  (0,1).  The  results  of  the  simulation  are  shown  on  Table 
6.1  for  the  gamma  distribution  and  on  Table  6.3  for  the  Weibull  distribution.  For  the  gamma 
distribution  we  chose  T  =  10  and  for  the  Weibull  we  chose  T  =  2. 

For  the  other  repair  model  we  will  use  the  supplemented-life  repair  model  with  0j  = 

Di  and  the  D/s  are  uniformly  distributed  on  (.8, 1).  We  will  refer  to  this  as  Model  II. 
In  the  context  of  our  interpretation  of  the  extended  life  repair  model,  restricting  the  T>j’s  to 
be  at  least  .8  restricts  the  increase  in  the  expected  remaining  life  to  be  at  most  25%  of  the 
original.  This  seems  to  be  a  reasonable  assumption,  hence,  the  choice  of  the  interval.  The 
results  of  the  simulation  are  shown  on  Table  6.2  for  the  gamma  distribution  and  on  Table 
6.4  for  the  Weibull  distribution.  We  use  the  same  value  for  T  as  in  the  case  of  Model  I. 

All  the  results  are  based  on  5,000  iterations  of  the  simulation.  To  generate  a  sample  for 
the  case  of  a  gamma  distribution  we  make  use  of  an  algorithm  by  Dagpunar(1978)  on  sam¬ 
pling  variates  from  a  truncated  gamma  distribution.  For  the  Weibull  case,  let  Zi,  Z2, . .  ■ 
be  a  random  sample  from  a  standard  exponential  distribution.  It  is  not  difficult  to  see  that 
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setting  Xj  =  (A“  +  ZjY^^  produces  the  desired  repair  process  for  a  Weibull  with  parameter 
a  under  Model  I  and  setting  Xj  =  (Zi  +  Z2  +  •  •  •  +  Zjf!^  produces  the  desired  repair  process 
for  Model  II.  To  generate  the  exponentials  we  used  the  function  REXP  given  by  Marsaglia 
and  Tsang  (1984). 

In  general,  the  simulation  seems  to  indicate  that  the  band  performs  well  in  both  cases 
with  as  low  as  a  sample  size  of  20  needed  to  get  within  2%  of  the  desired  confidence  level  in 
the  case  of  Model  I.  In  the  case  of  Model  II,  a  sample  size  of  50  is  sufficient  in  most  cases  to 
attain  the  same  accuracy.  The  discrepancy  is  due  mainly  to  the  fact  that,  in  view  of  the  way 
the  data  are  generated,  there  are  more  failure  times  per  sample  under  Model  I  than  there 
are  under  Model  II.  In  most  instances,  the  band  performed  better  as  F{T)  moved  further 
away  from  1  (i.e,  large  values  of  a  for  gamma  and  small  values  for  Weibull).  This  is  expected 
because  of  the  reliance  of  our  large  sample  results  on  the  assumption  that  F{T)  <  1.  In 
the  gamma  case,  this  could  also  be  attributed  to  the  low  efficiency  exhibited  by  the  data 
generating  process  when  obtaining  variates  close  to  the  tails  of  the  distribution.  Finally,  the 
results  seem  to  indicate  that  a  larger  sample  size  is  needed  to  attain  the  confidence  level  of 
99%  than  it  is  to  attain  either  a  90%  or  95%  confidence  level. 
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Table  6,1.  Coverage  Probabilities  of  100p%  Confidence  Bands 


for  Gamma(a)  under  Model  I 


n 

p  =  .90 

p  =  .95 

.. 

p  =  .99 

a  =  3.0 

P 

11 

o 

a  =  7.0 

a  =  3.0 

a  =  5.0 

O 

II 

a  =  3.0 

P 

11 

o 

10 

.8808 

.8990 

.9080 

.9242 

.9356 

.9520 

.9704 

.9766 

.9844 

20 

.8916 

.9080 

.9062 

.9352 

.9514 

.9522 

.9794 

.9848 

.9888 

30 

.8900 

.9056 

.9074 

.9366 

.9494 

.9516 

.9808 

.9852 

.9870 

60 

.8922 

.9000 

.9080 

.9394 

.9478 

.9552 

.9852 

.9878 

.9882 

100 

.8988 

.9058 

.9002 

.9464 

.9504 

.9474 

.9864 

.9874 

.9888 

200 

.8960 

.9020 

.9024 

.9462 

.9500 

.9512 

.9862 

.9906 

.9908 

Table  6.2.  Coverage  Probabilities  of  100p%  Confidence  Bands 
for  Gamma(Q')  under  Model  II 


n 

p  =  .90 

p  =  .95 

p  =  .99 

jj^HR 

a  =  5.0 

o 

11 

O 

CO 

II 

a  =  5.0 

o 

II 

a  =  3.0 

a  =  5.0 

a  =  7.0 

10 

.7406 

.7936 

.8478 

.7944 

.8448 

.8974 

.8732 

.9134 

.9506 

20 

.8146 

.8442 

.8838 

.8716 

.9024 

.9318 

.9364 

.9602 

.9730 

on 

.8722 

.8914 

.8988 

.9158 

.9368 

.9568 

.9714 

.9780 

on 

.8826 

.9036 

.9204 

,9338 

.9430 

.9750 

.9814 

.9854 

100 

.8878 

.8920 

.8934 

.9390 

.9444 

.9480 

.9850 

.9852 

.9898 

200 

.8986 

.8932 

! 

,9096  i 

.9410 

.9492 

.9500 

.9856 

.9896 

.9878 

>  ^ 


Table  6.3.  Coverage  Probabilities  of  100p%  Confidence  Bands 
for  Weibull(a)  under  Model  I 

n  p  =  .90  P  ~  ’95  p  =  *99 

a  =  1.0  a  =1.5  a  =  2.0  a  =  1.0  a  =  1.5  a  =  2.0  a  =  1.0  a  =  1.5  a  =  2.0 
10  .9044  .8934  .9018  .9452  .9400  .9406  .9792  .9774  .9770 

20  .8996  .8942  .9034  .9444  .9386  .9464  .9832  .9800  .9824 


^  I, 
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